Degenerate Bose liquid in a fluctuating gauge field 
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We study the effect of a strongly fluctuating gauge field on a degenerate Bose liquid, relevant to 
the charge degrees of freedom in doped Mott insulators. We find that the superfluidity is destroyed. 
The resulting metallic phase is studied using quantum Monte Carlo methods. Gauge fluctuations 
cause the boson world lines to retrace themselves. We examine how this world-line geometry affects 
the physical properties of the system. In particular, we find a transport relaxation rate of the order 
of 2k^T, consistent with the normal state of the cuprate superconductors. We also find that the 
density excitations of this model resemble that of the full t-J model. 
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Much attention has been focused on the normal state 
of the cuprate superconductors, which exhibit non- Fermi- 
liquid behavior over a wide range of temperatures. For 
instance, the resistivity is linear in temperature up to 
lOOOK in optimally doped La2_3;Sr3;Cu04. Optical con- 
ductivity measurements show a Drude peak with a scat- 
tering rate l/rtr of the order of 2A:bT. The Hall resistivity 
is suppressed relative to the classical value B/nec with 
a 1/T temperature dependence. The magnetoresistance 
is also anomalous in that it violates Kohler's rule [Q. 
Anderson Q emphasized the role of antiferromagnetic 
correlations in these doped Mott insulators, and intro- 
duced the idea of a resonating-valence-bond (RVB) state, 
which pointed to spin-charge separation at low energies. 
In a gauge-theoretic formulation |^,^, charge and spin 
are represented by bosonic holes ( "holons" ) and neutral 
spin-half fermions ( "spinous" ) interacting via a fluctuat- 
ing U(l) gauge field. 

In this Letter, we study the bosonic holons as a two- 
dimensional Bose liquid in the presence of a fluctuating 
perpendicular magnetic fleld. Physically, the interaction 
of a boson with the transverse gauge field, a, describes 
the motion of a vacancy in a spin background with fluc- 
tuating quantization axes [Q. The flux, $r, due to this 
gauge field represents the chirality of the spins around 
the plaquette r. We treat the gauge field in a quasistatic 
limit (which we justify below) where we consider an an- 
nealed average over static flux distributions which scat- 
ter the bosons elastically. (In contrast, Wheatley and 
coworkers emphasized the dissipative part of the scat- 
tering.) In the regime of strong gauge fluctuations, the 
boson world lines attempt to retrace themselves [||j7| , so 
that this system can be regarded as a bosonic analogue 
of the Brinkman-Rice problem [^,^. We find that su- 
perfluidity is destroyed even when there is significant ex- 
change among the bosons. This results in a degenerate 
Bose metal with interesting charge dynamics such as a 
Drude peak in the optical conductivity, consistent with 
the experimental scattering rate of 2k^T. 

We believe that our results for the boson model are 
relevant to the charge dynamics of the t-J model where 



electrons hop (with matrix element to) under a constraint 
of no double occupancy and with a nearest-neighbor spin 
exchange energy J. In the gauge theory for the uni- 
form RVB state in this model , gauge fluctuations at 
weak coupling are described by the correlator I?(k, w„) = 
(|a^(k,w„)p) = {xk"^ + 7|aj„|/fc)"\ where x is the or- 
bital susceptibility of the spinon fluid, 7 is the Landau 
damping coefficient, cj„ — 2?™//? and (3 ~ l/k-oT. We 
can see that these fluctuations are overdampcd with a 
relaxation time which diverges as 1/k^ . This justifies the 
quasistatic limit, i.e. D(k, a;,i) I?(k,a;„ = 0)(5„^0j 
a first approximation to the long-wavelength physics of 
this system. The flux distribution is then spatially un- 
correlated: ('I'r'I'r') = {b"^ / Px) ^r.r' where b is the lattice 
spacing. Note that, although ($^) decreases with tem- 
perature, the typical flux through a plaquette is of the 
order of a flux quantum (<f>o — he/ e) in the temperature 
range where the cuprates are normal: T / J > Td J ~ 0.1. 
A simple estimate, using a Fermi gas to calculate x, gives 
us ($2)1/2 _ (3T/7rJ)i/2$o > 0.25$o- A more detailed 
calculation gives ($2) 1/2 ^ o.5$o for T > 0.3 J [0. In 
the presence of such strong fluctuations, the behavior of 
the holons should be insensitive to the detailed value of 
(<I>2) and we will work with a large but constant ($2). 

To be precise, we study in this paper the model de- 
scribed by the imaginary-time action: S = S^+Sq where 

^B-/(E^l|^-^B(r)^dr (1) 
= ^,Y.D'\^MaO^M' = (2) 

k r ^ ' 

= -tJ2ie'^''blb, + h.c.) + ^J2n,(n,- 1) (3) 

with U/t^lona.nLxL lattice. The partition func- 
tion Z for N particles is obtained by summing over 
all periodic configurations of boson world lines {xc(r)} 
{0 < T < (3,a = 1, . . . ,N) and all flux distributions: 

a — 1 
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where is the action for the world hnes in the absence of 
a magnetic field. The path integration includes world-line 
configurations where the final boson positions, Xc(/3), are 
related to the initial ones, Xq(0), by a permutation of 
the particle labels, a. In other words, exchange may give 
rise to world-line loops containing more than one boson. 
Performing the annealed average over the Gaussian flux 
distribution, we obtain an effective action for the bosons 



alone: Scs = 



S2 with 



-y 



2/3 



L'(xq(t) - Xq/(t'))Xq-Xq/ drdr' (5) 
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where -D(x) = X^k^-^o ^C^' 0)^''^ " proportional to 
the Green's function of the lattice Laplacian. The k=0 
contribution is excluded because it represents a global 
change of gauge for a. 

The nonlocal current interaction ^2 can be better un- 
derstood in terms of the winding of the boson world 
lines around the plaquettes of the lattice. For in- 
stance, on an infinite plane, the phase factor in Eq. 

is exp(— 27ri ?«r$r/$o) where w-^ is the total 
winding number of all the world lines around plaque- 
tte r. Averaging over the flux distribution gives ^2 = 
(27r^($^)/$o) Sr ^r- To be more careful, we should 
restrict ourselves to distributions with zero total flux: 
$r = 0, which means that we should replace 



by Er 



(Xj. Wr) /L in the above formula. This 



winding-number sum has been termed the "Amperean 
area" of the world-line configuration We see that, 
in the presence of strong gauge fluctuations, the parti- 
tion function is dominated by paths with zero Amperean 
area. These are "retracing paths" where each traversal of 
a link on the lattice is retraced in the opposite direction 
at some point in time (Fig. |l|a), and so their contri- 
bution to Z is unaffected by the average over the gauge 
field. On the other hand, averaging causes destructive 
interference for non-retracing paths. 

It is in fact more efficient to enumerate S2 in terms of 
winding numbers. To do so under the periodic boundary 
conditions of a torus, we have to extend the definition of 
Wr. This can be done for world- line configurations which 
do not have a net wrapping around the torus in either 
spatial direction. A suitable definition (i.e. one which 
preserves Stokes' theorem, § a(x) • dx = w^^^: in the 
case of zero total fiux) is: Wj. = /[a°(x) — a^(x)] -dx, 
where a^ (x) is the vector potential at x due to a test fiux 
3" placed at plaquette r, and R is an arbitrary reference 
plaquette. Geometrically, this picks R to be on the "out- 
side" of any loop on the torus. 

We now turn to world-line configurations which have 
a net wrapping around the torus. These configurations 
are signals of superfluidity. Indeed, the superfluid den- 
sity rig (per site) is given by [^: rig = (W-^)/4/3f, 
where W = {Wx,Wy) are the total numbers of times 



the world lines wrap around the torus in the x and y di- 
rections respectively. In our problem, these paths pick 
up the random Aharonov-Bohm phases which thread 
the torus, and so should be suppressed after averag- 
ing over the gauge field. For instance, Eq. (ph gives 
52 = (1^,72/3) Efc.^o,.>=o^(k,0) - W^y{^^)L'^/^l for 
a straight-line path which wraps around the system Wy 
times in the j/-direction. In general, one can evaluate ^2 
for a world-line configuration with total wrapping num- 
bers W = {Wx,Wy) by decomposing it into a reference 
straight-line path with the same W and a set of world 
lines with no net wrapping (Fig. ^-c). Terms in Eq. (|^) 
involving the non- wrapping part only can then be evalu- 
ated as above, while terms involving the wrapping refer- 
ence path must be evaluated using Eq. d) directly. Thus, 
we see that any wrapping configuration would give a pos- 
itive contribution to 5*2 of order W^L^, and so should be 
strongly suppressed in the partition function. We there- 
fore argue that, for fixed ($^), superfiuidity is destroyed 
at all temperatures as L — > 00. 
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FIG. 1. (a) A retracing world line path projected onto the 
xy plane, (b) A wrapping path, (c) Decomposition of (b) into 
a reference path and a non-wrapping path. 

To establish the existence of this degenerate Bose 
metal and to characterize this phase, we have performed 
a path-integral quantum Monte Carlo simulation of 
this effective action with U /t — 4 at 25% boson filling 
with periodic boundary conditions. ^2 is evaluated using 
the prescription described above. Note that ^2 is real 
so that there is no sign problem in the path integration. 
Each Monte Carlo step involves the reconstruction of the 
world lines, {xq,(t)}, for all the particles (a = 1, . . . , N) 
in an interval of imaginary time. To ensure quantum 
exchange, we insist that each accepted configuration dif- 
fers from the previous one by a pair exchange. Since we 
are interested in the regime of strong flux fluctuations, 
we use a fine-grained discretization in the time direction 
(Ar < 0.1/i) so that deviations from a retracing path 
can be sampled correctly. A reasonably large number 
of imaginary-time points is also desirable for our ana- 
lytic continuation of the dynamical quantities of interest. 
These factors limit our simulations to T > O.li. 

We have measured the superfluid fraction &iT — t/Q 
for a range of flux variances and for system sizes up to 
7x7 lattices (Fig. ^). In the absence of gauge fields, the 
system has a Kosterlitz-Thouless transition at Tkt — i 
so that the superfluid fraction is nearly unity at T = t/Q. 
There is a threshold ($^)c beyond which the superfluid 
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density is exponentially small. Our arguments suggest 
that ($^)c vanishes as as L ^ oo (Fig. inset). 
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FIG. 2. Superfluid fraction as a function of flux variance 
($^) and system size L at /?t = 6 {tAr = 0.05). In- 
set: a scaling plot suggests that superfluidity vanishes at 
(-1'2)cLV$o ^ 8. 

Although our system has no superfluid response, it 
may still have diverging correlation lengths and diamag- 
netic response as T ^ 0, as in the classical analogue of 
this problem where the gauge field screens the vortex- 
antivortex binding potential and the Kosterlitz-Thouless 
phase is destroyed. Such divergences do not happen for 
strong gauge fluctuations. As already mentioned, boson 
paths are retracing in this limit and so, the system has no 
linear response to external magnetic fields. Indeed, the 
diamagnetic susceptibility is too small to measure in our 
simulation. This insensitivity to external fields is con- 
sistent with the observation that the Hall resistance is 
suppressed from its classical value, and that the magne- 
toresistance violates Kohler's law 

Another consequence of the retracing of paths is that 
the system might form dense aggregates. Feigelman et 
al. ]l3| showed that this phase separation might occur in 
the dilute limit. However, we have not found evidence 
for inhomogeneity or strong density fluctuations in our 
regime of moderate density and strong on-site repulsion. 

We now turn to the transport properties of this sys- 
tem. We use ("I>^) = 0.5$Q to ensure that we work in the 
regime of strong gauge fluctuations, where the world- line 
configuration have zero Amperean area. Using the gauge- 
invariant current j(T) = ^^[xQ(r + At) — XQ(T)]/Ar, we 
have measured the current correlation function to an ac- 
curacy of 0.3-0.9%. Through analytic continuation and 
the Kubo formula, the a.c. conductivity (y{Lo) is related 
to the current correlation function by 

-^o.M..(o))=y_^-^- (6) 

for T 7^ 0. We have performed this analytic continuation 
numerically using maximum entropy techniques p4[ . Us- 



ing the kinetic energy {K) , we can check, to an accuracy 
of 3%, the sum rule: a{uj)duj = -'k{K)IM?. For the 
lowest temperatures (T < 0.4i), we worked at fixed At 
and (3t/ L'^ with z — 2 to minimize finite-size effects [15|. 
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FIG. 3. Optical conductivity for 6x6 lattice at I3t = 
9,6,4,2,1,0.5. Inset: kinetic energy, proportional to spectral 
weight. Dashed line (BR) denotes Brinkman-Rice band edge. 
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FIG. 4. Transport. Resistivity (solid symbols) and scat- 
tering rate (hollow symbols) as a function of temperature. 
Scattering rate is deduced using the half-height of o{lo). 

We find that (j{uo) consists of a single Drude-like peak, 
which sharpens as the temperature is lowered (Fig. |^). 
The width of the peak gives us the transport scattering 
rate 1/Ttr. We find a temperature dependence consistent 
with: l/Ttr = CfceT with C = 1.8-2. (Fig. |). The resis- 
tivity p, given by the peak height, is consistent with a lin- 
ear temperature dependence of pe^/h = {l/2TrnQ)k-BT/t 
for T < 2t, (Fig. ^), where no is the boson number 
density. We estimate a statistical error of 5% for p by 
examining fiuctuations due to deviations in the current 
correlation function |16|. (There are also systematic er- 
rors due to the smoothing of structures.) There appears 
to be a systematic deviation from the linear-T behavior 
below T = 0.3t. This deviation is stronger for p than 
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for 1/rtr- The difference can be attributed to the T- 
dependence of the weight under the conductivity peak. 
Since we only have one peak which exhausts the sum rule, 
this spectral weight is proportional to —{K). As T de- 
creases, {K)/N drops below the band edge of —2^/3t for 
the single-particle problem Q and approaches remark- 
ably close to —4t per particle (Fig. || inset). 

When we examine the paths, we find that the world- 
line loops span increasing number of periods in imaginary 
time as the temperature is lowered, indicating strong par- 
ticle exchange below T = 0.5t. As T is reduced further, 
different loops begin to retrace each other. Thus, the 
transport properties shown in Fig. are clearly charac- 
teristic of a quantum Bose system. This should be con- 
trasted with the Brinkman-Rice result |^ for fceT ^ t, 
where the spectral weight decreases as (K) ~ T~^, and 
1 /rtr begins to saturate to a constant. This gives a linear- 
T resistivity which clearly has a different physical origin 
from that shown in Fig. ^. 

Our resistivity is in agreement, to within a factor of 2, 
with Jaklic and Prelovsek ||l^ who provided an approx- 
imate diagonalization of the t-J model on 4x4 lattices. 
They also find a Drude peak of width ^ 2k'QT. In ad- 
dition, they found a broad background which may be 
interpreted with a frequency-dependent scattering rate. 
This is absent in our boson model. This could be due to 
the absence of fermion degrees of freedom or our neglect 
of gauge-field relaxation. 
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FIG. 5. Dynamic structure factor. Solid lines: our results 
for 6x6 lattice at pt — 6 with t — 0.9to- Dashed lines: tj 
model with 4 holes in an 18-site cluster with to/ J = 2.5 |l8| . 

Finally, we examine the dynamic structure factor 
S{q,uj), related to the density correlation function by: 



K(r)n_q(0)> = / (e-™ + e-(^--)-)5(q,c.)dc. (7) 
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where nq(r) = J2a e^'^'^'''^'^\ This should be relevant to 
electron energy loss experiments. The loss of long-ranged 
phase coherence means that the structure factor does not 
sharp phonon peaks as in the superfluid phase. For fixed 



q, S'(q, has a broad peak as function of w. For q 
in the (tt, tt) direction, these peaks coincide remarkably 
with features found in the exact-diagonalization study of 
Eder et al. |l^] at a similar hole density (Fig. ||), after a 
moderate rescaling to t by a factor of 0.9. On the other 
hand, our results differ from those of Eder et al. in the 
(tt, 0) direction. These authors found two peaks in the 
structure factor whereas we find only one. 

In summary, we have studied a degenerate Bose system 
which is metallic due to elastic scattering with random 
gauge fields. We have demonstrated that many features 
of this model, such as the longitudinal transport time, 
indeed mimic the behavior of the full t-J model and the 
normal state of the cuprate superconductors. 
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